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Abstract 

Based on the geodesic equation in a static spherically symmetric 
metric we discuss the rotation curve and gravitational lensing. The 
rotation curve determines one function in the metric without assum- 
ing Einstein's equations. Then lensing is considered in the weak held 
approximation of general relativity. From the null geodesies we derive 
the lensing equation and corrections to it. 
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1 Introduction 



As long as the dark matter problem is open there is a non-zero probability 
that general relativity might not hold on the scale of galaxies [1-4]. Therefore 
a direct test on this scale is highly desired. It is the purpose of this paper to 
show how such a test is possible, if kinematical and lensing data of the galaxy 
are available. The idea is the following: the rotation curve determines part 
of the metric without assuming Einstein's field equations, only the geodesic 
equation is used. Then lensing can be calculated on the basis of the weak 
field approximation to general relativity and checked for consistency. In 
contrast to the usual way of analyzing the data no model for the galaxy 
must be constructed. This offers the possibility to test the basic physics. 

We treat lensing by means of the geodesic equation as well. By comput- 
ing the null geodesies we derive the lensing equation and we find corrections 
to it. Even if these corrections were not needed for the analysis of present 
day lensing data, they have to be under control for all eventualities. We 
only consider the static spherically symmetric case here in order to make 
the argument as simple as possible. 

2 Geodesic flow and rotation curve 

We consider a static spherically symmetric metric which we write in the 
form 

ds 2 = e v c 2 dt 2 - e x dr 2 - r 2 (cM 2 + sin 2 tidip 2 ) (2.1) 

where v and A are functions of r only. We take the coordinates x° = ct, 
x 1 = r, x 2 = x 3 = ip such that 

5oo = e", 5ii = -e A 

922 = ~r 2 , 533 = ~r 2 sin 2 (2.2) 

and zero otherwise. The components with upper indices are the inverse of 
this. The determinant comes out to be 

g = detg^ = -e u+x r 4 sin 2 (2.3) 

Let us recall the Christoffel symbols for the metric (2.1) from the ap- 
pendix of [3] 

r? = ^A ri = Ve-\ r} 1 = iv, v\ 2 = - re ~ x (2.4) 
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r 33 = -^- A sin 2 tf, r? 2 = i r§3 = - sin cos <?, r? 3 = i ri 3 = cottf 

and zero otherwise, the prime denotes the derivative with respect to r always. 
The geodesic equation is given by 

The origin of our reference frame is in the center of the galaxy. We consider 
geodesies in the plane 9 = ir/2, then we must solve the following three 
equations 

d 2 ct .dctdr 

+ u'—— = (2.6) 



ds 2 ds ds 

J2„ ,,/ 



(is 2 2 V ds ' 2\ds' \ds' 
d 2 cp 2 dr dip 



ds 2 ^ r ds ds ^ ^ 



Multiplying (2.6) by expi/ we find 

ds v as / 

so that 

,,<ici 

e — — = const. = a 
ds 



^ = ae~ v . (2.9) 
cis 

Next multiplying (2.8) by r 2 we get 

odip 

r — — = const. = J 

as 

where J is essentially the conserved angular momentum, hence 

dip J 
ds r 2 



(2.10) 



Finally, substituting (2.9) and (2.10) into (2.7) and multiplying by 2(exp A) 
xdr/ds we obtain 



ds 



\/dr\2 2 _„ J 2 
e ( — ) -a 2 e v + 



!/ *~ ' - " 0. (2.11) 



V ds ' r 2 
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Consequently, the square bracket is equal to another constant = b. Then 
the resulting differential equation can be written as 

(*)» =oV P*> + .-»(»-£). (2.12) 

To obtain the connection with the rotation curve which is an impor- 
tant astronomical observable, we remember the definition of the unitary 
4-velocity 

u - 



ds 

The term unitary indicates that u a has invariant length 1: 

U = 9 ^^s- = JaW = 1 - (2 ' 13) 
In our case u a is equal to 

\ ds ds ds ' 

Using (2.9) (2.10) and (2.12) we easily see that 

u 2 = -b=l (2.15) 

which by (2.13) fixes the constant of integration b = — 1. 

Clearly the last constant of integration a 2 must be related to the ge- 
ometry of the geodesies. To see this we consider the streamlines r = r(<p). 
Deviding (2.12) by J = r 2 dip/ds we obtain 

Introducing the variable 

«(„) = ;JL, (2.17) 



we write the equation in the form 

\dtp) ~ J 26 +e \J 2 



(™Y=^e-^) + e-H^- W 2 ). (2.18) 



To compare this equation with Newtonian dynamics we use the expansion 
of the metric for large r: 

e -(A+.) = 1 + (r- 2 ), e~ x = l-j + 0(r- 2 ). 
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Here 

2GM 

r s = —r- 2.19) 

cr 

is theSchwarzschild radius in case of a point mass. Then to order 1/r we 
have 

/dw\2 9 a 2 + b / b 9 \ , . 

(-) +v ? = ^— rs v,( 1 s-v?). (2.20) 

In Newtonian mechanics the bounded streamlines are ellipses 



w = - = -(l + ecos(p), (2.21) 
r p 

where p and e are parameter and eccentricity of the ellipse, p is connected 
with the non-relativistic angular momentum J by 

I 2 

— = GM. (2.22) 
p 

The Newtonian equation which corresponds to (2.20) now reads 

/dw\2 9 e 2 - 1 2 , 
(3-)+*= 5— + 2.23 

Comparing the coefficients in (2.20) and(2.23) we first find 



J L p 

By (2.22) and J = cJ this gives b = — 1 in agreement with (2.15). Secondly, 
from 

a 2 + b _ e 2 - 1 _ 
J 2 p 2 pa ' 

where a is the big half-axis of the ellipse, we obtain by (2.22) 



This shows that a 2 is connected with the big half-axis of the Kepler ellipse. 
The 3-velocity v which is measured by astronomers is defined as 

/dx 1 dx 2 dx 3 \ 

9 - (-*•-*•-*)■ (2 ' 25) 
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Using 



ds c 
dt a 



V 



(2.26) 



we can calculate 




(2.27) 



where u 2 is the spatial part in (2.13). Since our metric is diagonal it is 
simply given by 

_ 9 ^ dxi dxi n 
- 2 = E^^ = 1"« 2 ^- (^8) 

By (2.26) we now get the desired velocity squared 



As a check we determine the asymptotic behavior for large r 3> rs- 
Assuming circular motion (a = r) and using 1/a 2 = 1 + GM/cPr we find 



Since v = —r s /r we arrive at 




M is the total mass (normal plus dark) . This agrees with Newtonian dynam- 
ics (Kepler's third law). Summing up, the relation between observational 
quantities and theory is very direct. The rotation curve v(r) gives the metric 
function u(r) by solving the quadratic equation (2.29) 




(2.29) 




(2.30) 




(2.32) 



For velocities v -C c and r 3> r s this simplifies to 




(2.33) 
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3 Null geodesies and lensing 

In the case of null geodesies describing light rays the integration constant b 
in the geodesic equation (2.16) must be 

In the lensing problem one uses the weak field approximation to general 
relativity 

2 U(r) \ 2U(r) v , . 

where U(r) is the gravitational potential. The latter can be obtained from 
the rotation velocity according to (2.32) or (2.33). Expanding the square 
root in (2.33) for v 2 <C c 2 we get the very simple result 

U(r) = - V 4, (3.3) 
we see that the plus sign in (2.33) must be used. Introducing the quantity 

d=- (3.4) 
a 

in (3.1), the following first order equation remains to be solved 

The meaning of d becomes clear when we consider the trivial solution for 
17 = 0: 

d 

r = . 

sin ip 

It describes a straight line with distance d from the origin in polar coor- 
dinates (fig.l). After inversion the equation (3.5) can simply be solved by 
quadrature: 

dip ±1 



dr 



(3.6) 



The sign herein depends on the branch of the geodesic to be calculated. 
In case of a point-mass (Schwarzschild) lens we have 

U(r) - -™ (3.7) 
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and from (3.6) we get an elliptic integral for the polar angle (f(r): 

r 

^ " = ±d I ^-r>J r + ri)+rr a <P ' ^ 

where we have again used the Schwarzschild radius r s (2.19). To reduce this 
integral to Legendre's normal form we need the four zeros 01,02,03,04 of 
the quartic under the square root. We have the following four real roots 



a 1 



= d(Jl + j- £ -), o 2 =ed, (3.9) 



03 = 0, a 4 = d(-^l + j- |), (3-10) 

where 

e = 7 (3 - n) 

is a small parameter. It is convenient to expand everything in powers of e: 

2 2 
01 = d(l - J + ^-), a 2 = ed, o 3 = 0, a 4 = -d(l + J + ^-) (3.12) 

z o z o 

up to 0(e 3 ). The integral (3.8) is an incomplete elliptic integral of the first 
kind -F(<3?, k) where the parameter k is given by 

k*= ^- a ^- ai) =2e{l-e) (3.13) 
(ai - a 3 )(a 2 - a 4 ) 

(see [4], vol.11, p.310). 

As first application we compute the Einstein deflection angle and the 
correction to it. The origin of our coordinate system is at the mass M, 
polar axis goes from M to the observer (see fig.). We integrate (3.8) from 
the apex tq to infinity which gives us the deflection angle (p oo — ir/2. The 
apex is defined by the condition 

dip 

which gives ro = a± . Then we obtain 

¥>oo-|=/«*F($oo,fc), (3.14) 
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where the Jacobian fi is equal to 



P= J, \ S = -A 1 - £ o + t) ( 3 - 15 ) 

The argument follows from 

sin 2 = = l(l + h- 1^)= hi - cos(2<M) (3.16) 

(see [4], vol.11, p. 310). The elliptic integral can be expanded for small k as 
follows 



fc 2 1 



fc) = $ + —($- - sin2$) + 0(k A ) (3.17) 



(see [4], vol.11, p. 313). This finally gives 



< 3 ' 18 > 



The e = r s /<i is Einstein's result. 

Next we want to derive the lens equation. In this problem the observer 
is not at infinity but in a finite distance D d from the lens. The light source 
is at a distance D ds at the other side of the lens and an amount r] off the 
optical axis (see fig.); we use the same notation as in [5]. We have now 
to compute the null geodesies from the source at distance D ds through the 
apex r = ro,-d = ir/2 to the observer at distance D d . Then the polar angle 
j3 = -j^- of the source follows from 

7r-p = fid[F(^ ds ,k)+F^ d ,k)}. (3.19) 
Here the angle & d is given by 

. 2 « a 4 - a 2 D d - ai 

sin <Prf = = 

a 4 - cti D d - a 2 

= 1(1 _ cos(2d> d )) = i[i + l e -e(l- l ^e 2 ) + 9 2 £ ] (3.20) 

where Q = d/D d is the angle under which the observer sees the source. Then 
from (3.17) we obtain 

F ^^=4 +£ (2 + 8)-2-- + £ (8-8)- (3 - 21) 
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F(& ds , k) is given by the same formula with substituted by a = d/D ds . 
Now we find from (3.19) 

-P=-Q- a + 2e + e 2 (^-^)+0(e 3 ). (3.22) 
The lens equation is usually written in terms of the angles 

= P^1, 8 = -£-, a = 9^. (3.23) 
i-> s JJ d e JJ ds 

Then (3.22) gives the following lens equation with corrections 

a ~ D ds r s /3 TT\D ds / r s \2 

^ e - 2 ^D^-(-2-d-D:(iw) (3 - 24) 

for the point-mass lens. The first three terms are the leading order standard 
result ([6], p. 27). The third term can be written in terms of the deflection 
angle (3.18) as it is usually done. Comparing the corrections in (3.24) with 
those in (3.18) we find no direct correspondence. That means the lens equa- 
tion as usually written by means of the scaled deflection angle ([5]. p.21) is 
the leading approximation only. 

Now we turn to the formulation of the lens equation in an arbitrary 
spherically symmetric metric. Prom (3.6) we have the following integral for 
the polar angle of the null geodesies: 

r 

7T f dv 

^ -2= ±d J rVr'(l-u*)-(l+u)* ' ^ 

where we have introduced the dimensionless gravitational potential 

u(r) = ^U(r). (3.26) 

Here tq is the apex and it is important to note that only the potential values 
for r > ro contribute. For u = the trivial lens equation (3 = comes out, 
this follows from (3.24) for r s = 0. For u^O but |u| <C 1 the modification 
of the result comes from the neighborhood of the apex. It is therefore good 
enough to expand the potential u{r) in the vicinity of r = ro. For this 
purpose we use the beginning of the multipole expansion 



u(r ) = G0 + £L + £|. (3.27) 
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The constant term is necessary in view of the flat rotation curves; note that 
the potential has an absolute normalization in (3.2). It is unimportant that 
(3.27) breaks down for small r because we need r > tq only. 

With the three terms in (3.27) we get an elliptic integral of the first kind 
again: 

r 

Vir) ~ I = ±d f -£= = ±-^=F(*, k). (3.28) 



ro 



Here the quartic is given by 



G(r) = (l-cl)r i -2c c 1 r 3 -(d 2 +c d 2 +c 2 1 +2c C2)r 2 -(c 1 d 2 +2c 1 c 2 )r-c 2 d 2 -cl, 

(3.29) 

the Jacobian \x and the parameter k are the same as before (3.13) (3.15). 
The four zeros of G(r) are obtained by solving the two quadratic equations 



l + u(r) = 0, r 2 (l -u(r)) - d 2 = 0. 



This leads to 



ci d 2 + c 2 , cf 



2(1 -c ) V 1 ~ c o 4 ( 1 " c o) 2 



Cl cf c 2 

°2 = r + 



2(1 + c ) y 4(1 + co) 2 1 + co 



a3 -"2Tif^-V4Tif^-rf^ (3 - 30) 



ci I d 2 + c 2 c\ 

° 4 " 2(1 -c ) V !- c o 4(1 -c ) 2 ' 
Then the exact lens equation is contained in the analogous equation to (3.19) 

tt-P= ^ [F^ ds ,k) + F(^ d ,k)}. (3.31) 

The angles are given by the same formula (3.20) as before. The 

appropriate expansion of the lens equation (3.31) depends on the particular 
values Co,ci,C2 in (3.27). 

Regarding applications of our results one must replace the euclidean 
distances by angular diameter distances as usual. Galaxies with joint lensing 
and dynamical data can be found in the Sloan Lens ACS Survey (SLACS) 
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and its follow-up project [8]. Unfortunately, until today only one system 
SDSSJ 2321-097 has been analyzed in detail. This is an early- type elliptic 
galaxy which cannot be approximated by a spherically symmetric metric. 
So we must extend our model-independent analysis the the elliptical case or 
hope that the astronomers come up with a E0 lens galaxy. 
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